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Abstract. We report the first observations of the fourth harmonic {v^) in the 
azimuthal distribution of particles at RHIC. The measurement was done taking 
advantage of the large elliptic flow generated at RHIC. The integrated W4 is about 
a factor of 10 smaller than V2- For the sixth (ug) and eighth (ug) harmonics upper 
limits on the magnitudes are reported. 

Anisotropic flow, an anisotropy of the particle azimuthal distribution in momentum 
space with respect to the reaction plane, is a sensitive tool in the quest for the 
quark-gluon plasma and the understanding of bulk properties of the system created 
in ultrarelativistic nuclear collisions. It is commonly studied by measuring the Fourier 
harmonics of this distribution jj. Elliptic flow, V2, is well studied at RHIC and 
is thought to reflect conditions from the early time of the collision. Recently, Kolb j2] 
reported that the magnitude and even the sign of V4 are more sensitive than V2 to initial 
conditions in the hydrodynamic calculations. Besides one early measurement at the 
AGS IH] , reports of higher harmonics have not previously been published. Some of the 
present work has already appeared 

Experiment — The data come from the reaction Au + Au at a/snn = 200 GeV. 
The STAR detector main time projection chamber (TPC) was used in the analysis of 
two million events. The main TPC covered pseudorapidity (77) from -1.2 to 1.2 and the 
low transverse momentum (pt) cutoff was 0.15 GeV/c. In the present work all charged 
particles were analyzed, regardless of their particle type. The errors presented in the 
figures are statistical. 

Analysis — The difficulty is that the signal is small and the non-flow contribution 
to the two-particle azimuthal correlations can be larger than the correlations due to 
flow. To suppress the non-flow effects the current analysis uses the knowledge about 
the reaction plane derived from the large elliptic flow. One method for eliminating the 
non-flow contribution in a case when the reaction plane is known was proposed in pp. 
Results obtained with this method we designate by V4^{EP2}. The analysis for V4 was 
also done with three-particle cumulants jS] by measuring (cos(20a -|- 20{, — 40c))- 

Pt-dependence — The results as a function of pt are shown in Fig. ^ (left) for 
minimum bias collisions (0 — 80% centrality). Shown for are both the analysis relative 
to the second harmonic event plane, Vi{EP2}., and the three-particle cumulant, f4{3}. 
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Both methods determine the sign of V4 to be positive. As a function of pt, V4 rises more 
slowly from the origin than V2, but does flatten out at high pt like V2- The wefe) values 
are consistent with zero. Ollitrault has proposed [Hj for the higher harmonics that Vn 
might be proportional to V2^'^ if the distribution is a smooth, slowly varying function 
of cos(20). In order to test the applicability of this V2 scaling we have also plotted v| 
and V2 in the figure as dashed lines. The proportionality constant has been taken to be 
1.2 in order to fit the W4 data. The ratio, V4/V2, is shown in Fig. [T] (right) as a function 
of Pt. 




Figure 1. (left) The minimum bias values of V2, v^, and vq with respect to the second 
harmonic event plane as a function of pt for |7^|<1.2. The V2 values have been divided 
by a factor of two to fit on scale. Also shown are the three particle cumulant values 
(triangles) for V4 (v4{3}). The dashed curves are 1.2 • u| and 1.2 • iif- (right) The ratio 
V4/V2 is plotted against pt. The dashed line is at the value of 1.2. 



Parton coalescence — Assuming a simple parton coalescence model, for mesons one 
gets [7j 

v,/vl^l/A + l/2{vy{vlf). (1) 

Since experimentally this ratio is 1.2, v\ must be greater than zero. If one assumes that 
the hadronic scaling results from partonic w| scaling [S], then 

vl = {vlf (2) 

and 

v^/vl = 1/4 + 1/2 = 3/4. (3) 

But this is still less than 1.2. Therefore either vl is even greater than simple parton 
scaling would indicate, or the simple parton coalescence model is inadequate. 

Waist — Kolb j2] points out that for V2 > 10%, which occurs at high pt, and no other 
harmonics, the azimuthal distribution is not elliptic, but becomes "peanut" shaped. He 
calculates the amount of (which looks like a four-leaf clover) needed to eliminate this 
waist. Our values of as a function of pt are about a factor of two larger than needed 
to just eliminate the waist. 
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Centrality- dependence — The values of v^i^pt) for eight centrahty bins are shown in 
Fig. 121 (left). Integrating these values weighted with the yield gives Fig. |21 (right) which 
shows the centrality dependence of V2, W4, and vq with respect to the second harmonic 
event plane and also from three-particle cumulants {v^lS}). The Vq values are close 
to zero for all centralities. To again test the applicability of ^2'^^ scaling we have also 
plotted V2 and in the figure as dotted histograms. The proportionality constant 
has been taken to be 1.4 to approximately fit the V4 data. The larger constant here 
compared to that used in Fig. ^ is understood as coming from the use of the square 
of the average instead of the average of the square, and because the integrated values 
weighted by yield emphasize low pt, where the best factor is slightly larger. 




Figure 2. (left) v{pt) for the centrality bins (bottom to top) 5 to 10 % and 10, 20, 
30, 40, 50, 60, and 70 up to 80 %. (right) The pt- and rj- integrated values of V2, V4, 
and vq as a function of centrality. The V2 values have been divided by a factor of four 
to fit on scale. Also shown arc the three particle cumulant values for V4 (i'4{3}). The 
dotted histograms are 1.4 • v| and 1.4 • 

The Vn{EP2} values averaged over pt and 77 (|?7|<1.2), and also centrality (minimum 
bias, - 80%), are (in percent) V2 = 5.18 ± 0.01, = 0.44 ± 0.01, vq = 0.043 ± 0.037, 
and ^8 = —0.06 ± 0.14. Since Vq is essentially zero, we place a two sigma upper limit on 
vq of 0.1%. Also, vg is zero, but the error is larger because the sensitivity decreases as 
the harmonic order increases. 

Blast Wave fits — We have fitted the data with a modified Blast Wave model [Hj: 

p(0) = po(l + 2/2Cos(20) + 2/4Cos(40)) (4) 

^ Jl^d(f)cos{n(f))In{at)Ki{/3t){l + 2s2 cos(20) + 2s4 cos(40)) 

r.rf0/o(«t)i^i(A)(l + 2s2Cos(20) + 2s4Cos(40)) ' ^ ' 
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where J„ and Ki are modified Bessel functions, and at{(t>) = (pt/T) sinli(p(0)) and 
Pt{(l>) = (TTT't/T) cosh.{p{(())). In tliese equations, po is tlie transverse expansion rapidity 
(fo = tanli(po)) of the cyhndrical shell. The parameters /2 and are the harmonic 
amplitudes of the azimuthal variation of p, and S2 and S4 describe the spatial anisotropy 
of the source. 

The Blast Wave fits to V2 and V4 are shown in Fig. 01 (left) and expanded in Fig. El 
(right), showing the approximate agreement with the ratio. A temperature of 0.1 GeV 
was assumed giving the fit parameters po = 0.49, /2 = 1.4%, S2 = 9.1%, /4 = 0.0%, 
and S4 = 4.4%. It is interesting that in this large pt range the s values are considerably 
larger than the / values. 




p, (GeV/c) p, (GeV/c) 



Figure 3. (left) V2 and V4 as a function of pt witli tlic lines sliowing tlie Blast Wave 
fits, (right) The ratio V4/V2 as a function of pt with the line showing the ratio of the 
Blast Wave fits. 

Conclusions — We have measured V4 as a function of pt, and centrality. This is 
the first measurement of higher harmonics at RHIC. It is expected that these higher 
harmonics will be a sensitive test of the initial configuration of the system, since they 
provide a Fourier analysis of the shape in momentum space which can be related back 
to the initial shape in configuration space. 
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